ARKIV FOR FYSIK Band 15 nr 25 


Communicated 14 January 1959 by Oskar Kier and Ivar WALLER 


Transition matrix elements 
between states in a Coulomb field. II 


By Per O. M. Otsson 


Introduction 


In a previous paper on the same subject [1] a simple integral occurring in the 
calculation of transition probabilities between states in a Coulomb field was studied 
to some extent. It was felt that the ideas presented in the paper could be generalized 
so that essentially the same methods could be applied to a very general class of matrix 
elements in a systematic way. There was some difficulty in deciding upon the most 
promising course to follow. The final result presented in [1] is propably very nearly 
as simple as is possible and it seems to be advisable to maintain the structure of the 
previous formula as much as possible. It is worth noting that this formula can be 
obtained from one of the power series expansions for the Coulomb wave functions 
which was given in a previous paper [2]. The convergent expansion is then obtained 
after term-wise integration. This simple procedure, however, does not easily lead to 
the desired result if the ordinary expressions for the wave functions are used. These 
functions do not seem to be sufficiently manageable and therefore a representation of 
the wave function in terms of one of the confluent Appell functions has been in- 
troduced. This representation has increased symmetry and allows a great number of 
expansions and transformations. As is well-known, a very general class of matrix 
elements occurring in physics can be expressed directly in terms of one of the 
hypergeometric functions of several variables introduced by Lauricella. Examples of 
such matrix elements are the radiative transition matrix elements between states 
in a Coulomb field. Absence of the Coulomb interaction is of course obtained as 
a special case provided that the interaction, which is then non-electromagnetic, can 
be expressed in suitable terms. 

The introduction of the Appell function representation allows, as mentioned above, 
a great number of expansions and transformations. Such expansions and transforma- 
tions lead to a great number of analytic continuations of the Lauricella function, 
in which the Lauricella function is, in general, expressed in terms of other hypergeo- 
metric functions. These expressions are as a rule complicated and it is also in many 
cases very difficult to decide whether the expansions are convergent or not. 

It was strongly felt that, even in the more general cases when the matrix elements 
lead to a Lauricella function instead of an ordinary hypergeometric function as 
was the case in [1], the matrix elements could be expressed as a sum over terms con- 
nected by a recurrence relation in a similar way to that in the previous paper and that 
such an expansion of the matrix element should be convergent without modification 
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for all values of the parameters which are of physical interest. For those values of the 
parameters involved for which the expansion given in [1] is too slowly convergent for 
practical purposes, the expansion can be easily modified so that the rapidity of the 
convergence is increased without changing the structure of the formula. This is of 
importance if an electronic computer is used. It is desirable that this property is 
maintained also for the more general types of matrix elements which may occur. 
This means that analytic continuations of the Lauricella functions should be given 
which are not only convergent, but also sufficiently rapidly convergent for practical 
purposes, for the whole domain of values of the arguments of the function which 
are of physical interest. Here, the selection of results of physical interest is of 
course only a tentative one. 

Although the resulting formula for the matrix element obtained in [1] was not 
obtained from term-by-term integration after performing a series expansion of the 
integrand, it could, as was mentioned above, have been obtained in that way. 
When using such a procedure the convergence of the result depends on the choice 
of a proper expansion of the integrand. For the more general matrix elements that 
will be considered in this paper, the known expansions of the Lauricella function are 
not convergent or are very complicated whereas the expansion of the integrand is 
greatly simplified by the Appell function representation of the wave functions. The 
problem of finding useful expansions of the Lauricella function is then conveniently 
replaced by the problem of finding a proper expansion of a product of confluent 
Appell functions. 

Since the simple power series expansion led to a satisfactory result for the integral 
studied in [1], it is reasonable to try the same method also for the more general types 
of matrix elements which will be discussed in this paper. In the following, it will be 
shown that the Lauricella function, and thus all matrix elements that may be ex- 
pressed in terms of it, can be expanded in a series that is generally convergent and 
whose terms are connected by a recurrence relation. The number of coupled fields 
is not important but the recurrence relation connecting the terms of the expansion 
will become more complex and involve an increasing number of terms with 
an increasing number of coupled fields. 


An expansion for products of fields and their matrix elements 


The matrix elements that will be studied can be written 


- N 
I= eden TI ®, (0, 653 e135 a7, yr) dr, (1) 
rj i=1 
where ®, is a confluent Appell function defined by the convergent expansion 
= (b)m (B')nv” y” 


0 (C)minm! n! ° 


®, (b, b'; c; @, y)= (2) 


It seems to be sufficient to connect the pairs of parameters b; and b; by the relation 

b, +b; =c; in which case the Appell function is reduced to more elementary func- 
tions [3]: 

@,(b,b';¢;,y) = 6" PF, (b;c5-y). (3) 
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From (3) follows the property 
e’D, (b,b';b + B32, y) =O, (b,0'5b +052 +A,y +A), (4) 


which means that we may put 4“ = 0 in (1) without loss of generality. We thus have 
a certain freedom of choice of the arguments 2,r and y;r in (1) and, in fact, each of 
the many possible expansions of the ®,-functions with the restriction b; + bj =, will 
contain an arbitrary parameter having an influence on the convergence of the expan- 
sion. The integral (1), containing the reduced ®,-functions (3), has been given by A. 
Erdélyi [4]: 


ve N 
ee [ete Tle” Fi (bs; 3 (u—y) "dr = 
t=] 
0 


T(a 
7 (a) area Oe Dette can Ons Oy cou Cns. Gig Sas nee Ch) (5) 
(oe p ») 
i=1 
where jee ie te UE 
fei > vi 
Tat 


and F’, is one of the Lauricella functions of N variables. It is defined in the neighbour- 
hood of the origin by the expansion 


ED eg... Ew, (6) 


It is not easy to find analytic continuations of this function beyond the domain 
of convergence which are of satisfactory simplicity. We therefore consider (5) as an 
integral representation of the /',-function from which continuations will be derived. 
From the previous paper [1], where the special case VN =1 was studied, we know 
that a satisfactory result can be obtained by expanding the ®,-function in powers of 
r. We shall therefore expand the N-fold product of Appell functions in the integrand 
of (1) in powers of r. In order to do this, we shall use an integral representation of 
the ®,-function [2] 


®, (b, b'; ¢; x, i 


ee hey 


“POPC ae = le Tu+y(l—u) ue- 11)? |F [e—b—b'; c—b; y(w—l)] du. 


ds 


If we use throughout } + b’ =c the first parameter of the , F,-function is zero and 
thus : 
[T(b+0’) | +y@Q—u),,b-1 b/-1 
oF ‘¥ a ACA ru u te p 
®, (b, b'; b+b'; x, y) Pore), e w*(l—u) du (7) 
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The product of two such functions can then be written 


D, (b,, bi; b, +b; 17, yy 1) Dz (bo, bs ; b, + b3; LyX, Yr) = 


1 
— Lr (b, + bi) (r7,u+ y¥,Q—)jr He , fee ni 
- POI CD , ¢ ®, (b,, be ; Ae Maps sce Yot) U (1— ) du 


1 
T (6, +03) " es ; 
“PT 6) Jo *(1—u)?* @, [b,, 62; by +02; (w+ xu + y, (1— u)) 
(Yo taut y (l—u))r]du. 


We have used here the transformation (4). The ®,-function can be expanded 


2[b.; bz; by + b3; (x2 +x, uty, (l—wu))7r, (Yat uuty, (1—u))r] 


n = 
(y+ Yo)” 2-24] ” oe 
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— }-4-+y 
“+ Ye 


n=0 


This expansion is derived in [2]. Inserted into the previous integral it yields 


®, (b,, b:3 b, +013 27, Yy 7) Dy (bg, 623 by +023 XT, Yo7) 
1 
fo ; — n 
1 (b 1 + bi) = lta + vale jet a—a (1-4 * 4) - 
Wt 


T(b,) (81) no J 
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But the F',-function of two variables, which is usually called an Appell function 


ee 
and denoted F,, has the integral representation [5] 


a ee P(c)T(c’) 
F, (a, b,b > €,€ ; % =e) —b) 
ae 
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00 
Using an integral representation of the ordinary hypergeometric functions 


x 
E . 
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the integration over wu can be performed, yielding 


F, (a, b, b'; c, c'; x, y)= 


5: 3882S eg ee b’-1 — 4)\¢’—0’-1 7 —a bf te x 
rere =n |" (L—») (U0 0)) arg Ley (2.3; 0: = -\ae, (10) 


Comparing this integral representation of the F,-function with the coefficients 
of the expansion (9), it is easily seen that the latter can be expressed in terms of the 
F,-functions: 


®, (b,, 613 b, +b1; 27, y, 7) Dy (bg, O23 6, +02; %Q7, Yor) = 


= [Yr + Yo) 7)” ( Y, — &. , a=) 
F b,, bg; 6, + bi, 6, +62; +—*, 4 —*}. (11 
3, n! AIR a 5 ae Perr re Be 


Since the first parameter in the /’,-function is a non-positive integer, the F',-func- 
tion is a polynomial in its arguments. It is further obvious that the expansion is 
convergent for all finite values of r since it must have the same radius of convergence 
as the ®,-functions. 

It is now not difficult to see that the product of an arbitrary number of confluent 
Appell functions can be expanded in an analogous way: 


N 
T] ®2 (0: 0:5 b:+0;; air, yr) = 
i-1 


a AD ETA YX Yn — Ly 
_ INO tcc (One Craton CH renee ; 2 
»3 A F, Oi) cian Oye Ce Cy Dy, Dy, (12) 


where c;= 5, + bj. 


The restriction 5; +b; =c, is not necessary but is convenient. Similar expansions 
exist with the /,-functions replaced by more general hypergeometric functions 
but such functions have not been studied. Since the task of deriving the necessary 
properties of the functions, which is E Nece laborious, has to precede their use, 
we will only consider the case b, + b; 

The proof of (12) proceeds sintlacly ‘the proof for the product of two Appell 
functions. The integral representation (7) is introduced for all the ®,-functions but 
one and this latter function absorbs the exponential factors Seourting to (4) before 
it is expanded. The coefficients will then be N — 1-fold integrals, in which one ,F- 
function will occur in the integrand. These integrals are integral representations of 
the Lauricella F',-function of NV variables in virtue of the integral representation [6] 


Bg (By a Dey nc Oud Cys Cyne Cys Ly, Vqy.0-- Un) = 
= LT (e,) T (eg) 1 (ey) ee 
I (6,)T (6,) ... T (6y) T (c, — 5,) T (c,— 5g) ... P (ey — by) 
Lu 1 
x{ | aS [ atetuktaye day rt deayent Setl—apy ys NT x 
00 é 
x (1 — Uy — Up %y— 11° — Uy ty) “du, ...duy (13) 
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which, with the aid of the integral representation of the ordinary hypergeometric 
functions, can be written 


BAG 04; Obs a0 On 3s, CiscCasce nw Cry sattrs clone ony) = 
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Inserting the expansion derived into the integral (1) we obtain 


ie WOO! 
% Byil tue Bay; Dayads clbyiibyt Bebace Dement ba elon a ee ee 
2 24 


ioe) 
x [ er retetar= 
0 
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#4 uw n=0 n! s 
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When it is convergent, this series gives the integral J. The convergence depends 
mainly on the values of the parameters 2;,, y,; and «4. We now observe the possibility 
of adding arbitrary quantities to #,, y; and according to (4). This invariance of the 
integral J does not manifest itself in the F',-function but affects the convergence 
of (14) and results in an arbitrariness in the expansion of the integral J. Using (5) 
we may express this circumstance as an expansion of the Lauricella function in which 
one free parameter occurs: 


oe, ame 82y\ | 
1+s’l+s’ (1+8 


F, (« Dy, Oa, -:. Ons Cy, Cgs ao Ons 


oo a)n ‘ 
=(1+s)* 2, On (— 6) Fy %, Og, .2. By 5c, .'.. Ons Sas Sep woo Re OD 
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New parameters have been introduced in order to obtain a result that is more 
easily surveyed. 

A great number of expansions of the Lauricella function can be derived starting 
from (12). As an example we derive an expansion containing products of Lauricella 
functions with fewer variables. 

We have 


N 
TI ®, (6, 63; 6; 4+6;; ar, yr) 
re! 


N 


2 (Bi, bi; b+ bi; UT, yr) I] OD, (b;, 613 b:+0;; ar, mT) = 
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ll 
ey 
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Using this expansion instead of (12) we obtain for the integral (1): 
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From this expansion we obtain an expansion of the Lauricella function in a similar 
way to that in which (15) was derived 


SZ SZ 8Z tz tz 
1 IE ea ocean = Be ee ey eaeldal v_)— 
af Loa g ty wae "*L+eet? Leeeh Leet Al Peer al ee 
ae (2) m+n m n 
=(1+s+2) rE, Pamela (—+t) (1 Dae OCs maeaCiosteys ose h ee 
med s H ; 
XFy(—7, Div age by; Chetty +++ ONS 2K44> «++ Zn). (18) 


The summation over n can be performed according to (15) and the result can be 
expressed in the following way: 


eae ers le 
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v4 
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=(1+ 8)*> mI (18) FA = 1520), Oy.) -<- (0g) C1.) Caicos Ces eee ey 
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One parameter is lost during the summation. 
A special case of (19) is obtained if we put z, =z. =2, = ++ =% =0. For this case, 
we obtain an expansion analogous to (15) 


Zz z z 
FE, (« Dy, Ome sen On suGinCencses (ONS re eee sf i) 
= (1 a = (@)n ni b,.b b ‘ = 9 
=(1+ 8) 2, at (=8)" I, (a@-+'n, b,, 0, ... Byi Cy, Coy vee ONS Sys 29s ookew le (20) 
n= : 


For N =1, the formulas (15) and (20) are expansions of the ordinary hypergeo- 
metric functions. These special cases are found in Erdélyi, Magnus, Oberhettinger 
and Tricomi, Higher Transcendental Functions [7] and can be derived from certain 
bilinear relations between ordinary hypergeometric functions which have been given 
by Meixner [8]. Some of the relations given by Meixner follow from the expansions 
given above with N = 2 using an expansion of the F,-function in terms of products 
of two ,F,-functions, that has been given by Burchnall and Chaundy [9]. Erdélyi 
has also given a number of expansions of the ®,-function [10] from which a great 
number of other expansions of the Lauricella function may be derived. The expansions 
given here have been chosen because of their simplicity and because they are well 
suited for computations. By a proper choice of the free parameters in the expansions, 
it seems to be possible to obtain analytic continuations of the Lauricella function 
in a great number of cases of physical interest. The existence of recurrence relations 
for the F',-functions makes a numerical computation very convenient and the com- 
puting time will be proportional to the number of terms which are used to approxi- 
mate the /',-function. 
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Applications 


As a first application of the results derived in the previous section, we will consider 
the simple but instructive case when only one ®,-function occurs in the integrand: 


I= Je*r® *O,(y+in, y—in; 2y;ipr, —ipr)dr. (21) 


SS es Bs 


This integral is a special case of an integral that was considered in a previous 
paper [1] which contained one ®,-function and for which a convergent expansion 
was derived without the restriction b + b’ = ¢ that reduces the ®,-function to ordinary 
confluent hypergeometric functions. Because of the simplicity of the integral, it is 
well suited as a first example to which we can apply the previously derived 
expansions. It is also instructive to compare the following derivation with the result 
obtained in the previous paper. 

Using the results of the previous section, we immediately obtain 


T (a) iP J 
see -Qy: a 
Ii aap ehs (a 96: Wea née 
IT (a) = (B)n 2 , 
<page lt eS SR (aR (mms ytins rie) (2) 
where ee ee 
l+s 1l+itp 


We will assume that all the parameters occurring in (21) are real positive quantities 
except 7 which may be negative or zero. From the definition (2) of the ®,-function, 
it follows that it is a real function: 


®,(y +in,y — 1; 2y;2ipr, —ipr) =D4(y —in,y +1; 2y; —ipr,ipr). 
This means that we also have 


Tai 
 (l-ip 


—21p 
a (ay Uys ee =, 


from which the singular nature of the integral for 7p = + 1 is clearly evident. This 
can also be seen from the asymptotic behaviour of the integrand of the integral (21). 
It is then natural to choose the factor in front of the expansion (22) so that it takes 
account of the singular behaviour of the integral as much as possible. A good choice 
should be 


1 
l+8=7 oe 
We then obtain 
Oe act iD BA wey = 
Perey “ 1+ip’ 


297 


P. O. M. OLSSON, Transition Matrix Elements 


and thus 


P (a) 3 Gr ( —tp 


: 2 
_ Sais Daye 9 
(1+ p*)* nao n! soa 2F,( nm ytin; 2y; ), (23) 


lip 


This expansion converges for all real p-values. In order to prove this statement, it 
is sufficient to show that the ,/',-function is absolutely bounded: 


¢ 


\\<4 all n. 


This follows from the integral representation 
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2 2u_\" 
- ) const fw ay (1 — s du, 
l+ip 
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9 
since 1 — 


The above choice is not the only one possible for s, but it has the advantage of 
providing an expansion for which the convergence is easily proved as well as a 
natural method for obtaining the value of the parameter s. 

The expansion (23) can be written 


pu Ole 2 
T= BP eo ee 


where the c, are connected by a recurrence relation 


atn Be 
Cnt (2 +n) (n+1)1+p 


3 {2[p(n+y)—y]en—p(atn—1)¢n-y}, (25) 


which can be easily derived from the recurrence relations for the ,/,-functions. 
The ®,-function in (21) reduces according to (3) to confluent hypergeometric 
functions: 


—ipr 


O,(y+in, y—tn; 2y;ipr, —tpr) =e Fy (y+in; 2y; 2ipr). 

Expansions like (24) for the integral (21) can be obtained without this restriction 
but, when the integrand contains more than one ®,-function, we shall consider func- 
tions with this restriction only, since a study of integrals over products of general 
@,-functions seems to require the introduction of more general hypergeometric 
functions than the Lauricella and Appell functions. 

As an example of a ®,-function with general parameters, we give an expansion in 
,/’,-functions: 


Que 


®, (a,b—a;c; x+y, y)= > 
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The function can also be expanded in cylindrical functions [2]. 

The integrand of (21) can thus contain arbitrary confluent hypergeometric functions 
as well as certain expansions in such functions and cylindrical functions. It is not 
our purpose to investigate to what extent the expansion (22) converges in such 
general cases. It will be considered sufficient that a convergent expansion can be 
derived in cases where the parameters and arguments have values that are represent- 
ative for a great number of physical applications. 

It is interesting to observe that also for large p-values the formula (22) allows 
an expansion that seems to be sufficiently rapidly convergent for practical purposes. 
This indicates a rather general applicability of the method. 

In the limit p— co, or if the exponential factor e~” is not present in the integrand, 
we have to compute the ,F,-function of argument 2: 


: : a ~ (rn : ; 
oF, (a, y+in; 2y; 2)=(1+s) Pe nl (8) aFi(—m y tin; 2y; 2) (26) 
with estas 8) 
1+ 
We choose z=10, @ real, 
2 
then Aang ers 


If |@|<1 it can be shown that [11]: 


I 
Fi (—n, y+in; 2y;%0)|=0 = 
|2F1(—m, y+in; 2y; %9)| laser 


which proves the convergence of (26). 
As a second application we shall consider the integral: 


~ 2 
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where 


Both the ®,-functions in (27) reduce to ordinary confluent hypergeometric func- 
tions according to (3). From the asymptotic behaviour of the , /',-function for a large 
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argument, the singular nature of the integral for i(p, + p,) = + 1 can be easily seen. 
We therefore put 


l+s 1 1 
: : ; - : =— 2, 
1+i(p,+ py) [1+4(p,+P,)] [1480p =p, -8(p, 4011-4 (7,- 9) 1th 


(p, — Po)" = 4 (p, + Ms) [1 + (p, — Pe)*] 
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Thus |s|<1 for real, finite values of p, and py. 
20 py 
Further ON a en 
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For these values of z, and 2, F,-function in the expansion (28) is absolutely 
bounded: 


| F',( = 7, D1,0g; by Orso 20252], 24) |= A, el 


This can be shown from the integral representation of the F',-function given in 
the previous section, 
F,(—7n, by, 6,3; 6, + bi; bg +623 %, 29) 
1 1 - , , 
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h—1(p, + po) 
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since | 1 — uw, 2, — Uy 2,| 


for 0 <u, <1, 0<u,<1 and p, and p, real and positive which makes h real as a 
consequence of (29). For the sake of simplicity we assume that the parameters a, b; 
and 6; have such values that the integral representation used for the F',-function 
actually has a meaning. 

The procedure outlined above can obviously be applied to integrals containing arbi- 
trarily many ®,-functions, but the rapidity of the convergence seems to decrease with 
an increasing number of ®,-functions in the product for the given choice of s. In addi- 
tion, the recurrence relations for the F',-functions in the coefficients of the expansions 
will become more complicated and the number of functions connected by the recur- 
rence relation will increase. For the integral (27), where /’,-functions occur in the 
terms of its expansion, no more than five terms are connected by a recurrence rela- 
tion since it is known that five #’,-functions, where the parameters differ by integers, 
are connected by a linear relation in which the coefficients are polynomials in the 
arguments. 

The recurrence relation for the /,-function with integer steps in the first para- 
meter does not seem to have been derived. It is, however, possible to avoid the diffi- 
culty of obtaining recursion relations for the /',-functions by using expansions into 
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products of F',-functions with fewer variables, such as (19). As an example we may 
consider the expansion (28). Since it is absolutely convergent it follows that the use 
of the expansion (19) with the same value for the parameter s will also be convergent, 
since the series may be summed in an arbitrary manner. Instead of using the recurrence 
relation for the F,-function we can then use the simpler recurrence relations for 
the ,F,-function, with the unimportant complication of using it twice for each term. 

The general nature of the integrals investigated in this paper makes it practically 
impossible to state to what extent the given expansions can be used to compute the 
integrals. From the two applications above, it is obvious, however, that simple and 
rapid methods for computing integrals of the type studied can be obtained in a great 
number of cases of physical interest. 
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